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Abstract 

For flows of a lossless gas, the stagnation enthalpy obeys a lin- 
ear convected wave equation with coefficients which depend on the 
flow variables. This equation is self-adjoint and one has a reciprocity 
relation between source and observer. It fulfills for subsonic flow a 
quadratic conservation equation implying stability. It is taken as ba- 
sis for an acoustic analogy and is applied to the sound generation by 
the collision of a convected vortex and a rigid cylinder. 



1 Introduction 

The concept of an acoustic analogy was introduced by Lighthill pp. This idea 
offered him the possibility to derive important results on the generation of 
aerodynamic noise without relying on expansions or perturbation schemes. 
This is especially important, as the acoustic energy generated by an unsteady 
flow field represents usually only a minute part of the energy flux occurring 
in the flow and errors which are small compared to the flow quantities may be 
very large if one compares them with the sound quantities. The conditions 
under which perturbation schemes based on small sound level or on small flow 
Mach number may be used, have been clarified in later researches [21 El S] , 
although there still remain open questions. 

Lighthill derived without any approximations an exactly valid equation, 
which admits an acoustic interpretation. He showed that every compressible 
flow fulfills an inhomogeneous wave equation with a quadrupole type source 
distribution. As the inhomogeneous wave equation describes the generation 
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and propagation of sound in an ideal acoustic medium at rest, he had related 
aerodynamic noise to sound waves in the ideal acoustic medium. This relation 
is exactly true, but it relates the effects of convection and refraction by a 
steady basic flow to a source distribution in an acoustic medium without 
flow and this is difficult to visualize. One is therefore inclined to look for a 
relation to an ideal acoustic medium in motion. 

Here we will identify the operator which describes the propagation of 
sound in irrotational, isentropic flow and we will base the analogy on that 
operator. This of course does not mean, that vortices or entropy inhomo- 
geneities are excluded. It is very similar to the situation found in Lighthill's 
analogy. There, flows were exluded from the medium and they occurred as 
sources. Here, vortices and entropy inhomogeneities are excluded and they 
occur as sources. The equation is therefore very well suited to study the 
sound radiation from vortices which are convected in an irrotational flow. 
The acoustic variable which we use is the stagnation enthalpy. This vari- 
able has been used with aeroacoustic applications in mind before [5j [6] . The 
equations considered before differ however from that derived here, but com- 
pare [TJ. The operator which replaces the wave operator here is a self-adjoint 
one, even if the flow field is completely arbitrary. It is not necessary that 
it is irrotational or that it fulfills the basic eqations of fluid dynamics. The 
self-adjointness leads to a reciprocity principle which is valid for an exchange 
of source and observer. It also implies the existence of a variational principle 
from which the basic equation of the analogy can be derived. This principle is 
the simplest generalization of the variational principle for the wave equation, 
partial derivatives with respect to time are replaced by material derivatives, 
which seems natural if one requires Galilei-invariance. From the variational 
principle one can then conclude the validity of an energy theorem. The main 
drawback is of course the lack of a simple explicit expression for the Green's 
function. Notice however, that for certain important situations a Green's 
function to first order has been determined in [8|. An alternative would of 
course be a numerical solution. To check its feasibility, a general purpose 
PDE-solver for a PC has been applied to calculate the sound generated by a 
two-dimensional vortex convected along a circular cylinder. 

2 The acoustic analogy 
2.1 Preliminaries 

Lighthill based his theory of aerodynamic sound generation on an eqaution, 
which he obtained by cross-differentiation from the Euler equations and from 
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the continuity equation 

dpv 

namely 



V ■ pv v + Vp 
dp ^ 



cr at" 1 

Here p denotes the density, p the pressure, a the speed of sound, and v the 
particle velocity. This equation is valid for isentropic flow if losses and tem- 
poral variations of the speed of sound can be ignored. Temporal averages can 
be subtracted, therefore one may assume that only fluctuating quantities are 
contained in eq.fll]). No linearization has been performed in the derivation. 
Therefore eq.fll]) is valid under very general conditions. Often, especially for 
low Mach number flows, one neglects the acoustic contributions in the double 
divergence in eq.fll]) and considers this term as a known source term for the 
sound generation. Then one describes the sound field as a wave obeying the 
wave equation, i.e. as a sound field propagating in a non-moving medium. 
Convection and refraction effects are then neglected. Eq.fll]) is however gen- 
erally true and these effects are in principle contained in eq.fll]). Attempts 
to extract them have been made by separating the velocity into a mean and 
fluctuating part and thereby to obtain these effects. Although it should be 
possible to describe mean flow effects in this way, it has been felt that a more 
apropriate description should be obtained through a modification of eq.fll]). 
Instead of the wave operator, a "convected wave operator" which describes 
the propagation of sound in a moving medium seems more apropriate. One 
then rewrites all fluid quantities clS db superposition 

V -)> v + v,p -)> po +p, ■ • • 

of " nonacoustic" and "acoustic" variables and rewrites the basic equations 
in terms of these variables. If one neglects contributions which are nonlinear 
in the non-subscribed variables, one obtains the acoustic equations. It is well 
known, that one can derive from this system one " convected wave equation" 
for unidirectional (in x-direction) shear flow or for potential flows, namely 

1 D „ 2 _ 1 D 3 p ,„ , „ dp D d „ 

£shea ^ = ^^ v ' GoVp -^^- 2(VMo) ' v ^ = Wlth m = di + ^ 

for unidirectional and with an acoustic potential 

, p D<f) d p D(f) 

C pot <P = V ■ ( PO V0 - ^— v ) - = 
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for irrotational flows. For aerodynamic noise one should not neglect the 
nonlinear terms. It is however possible, very similar to Lighthill's approach, 
to collect the linear and nonlinear terms and to derive an inhomogeneous wave 
equation. It is even possible - and it has been done e. g. by Tarn and Auriault 
[9] - to follow this approach in the full continuity and Euler equations and 
derive an inhomogeneous linear system of equations. Then one could use for 
v , po, etc. the temporal mean values. The equations which are obtained 
by this method are much more complicated than the wave equation and are 
usually solvable only numerically. Furthermore they are rather different from 
the equations usually studied in mathematical physics and little is known 
about existence and uniqueness of solutions. It is however known, that the 
linear parts of the equations agree with the stability equations and as many 
flows are unstable, one has to expect, that the equations will be unstable. 
Of course, exponentially growing solutions are physically excluded - at least 
for longer times - and therefore they do not occur in the correct solutions. 
Small errors will however produce these instabilities and provisions have to 
be made to limit their growth. How these provisions influence the sound 
obtained from these calculations is difficult to assess. 

Here we follow a different strategy which is related to the propagation of 
sound waves in potential flows, but differs from the above described method 
significantly. We do not separate the velocity into an irrotational part and a 
remainder but follow a path used by Howe [5]. He observed that Bernoulli's 
equation states that the stagnation enthalpy differs from the potential only 
by a sign and by a time derivative. For irrotational sound waves one could 
therefore use the stagnation enthalpy instead of the potential. The stagnation 
enthalpy is however defined also for rotational flows and could be used as a 
variable of an acoustic analogy for arbitrary flows. This is what we will 
do. We will derive a convected wave equation for the stagnation enthalpy. 
We will discuss its main properties and show that it possesses many of the 
formal properties of the ordinary wave operator. An important example is 
an energy conservation theorem with an energy density which is positive 
for subsonic flows. This excludes instabilities, the influence of small errors 
remains small. The provisions which are necessary in many other analogies 
to limit the troublesome growth of instabilities are not necessary here. 

2.2 Basic Relations 

To obtain the equations for the acoustic analogy including convection effects, 
one starts from the Euler equation for compressible flow. Crocco's form of 



4 



these equations reads 



<9v 

+ V5 = -L, L = wxv-TVs, w = curlv. (2) 

at 

B denotes the total enthalpy B = /i+|v 2 with the enthalpy h and the velocity 
v. T is the temperature und s the entropy. From the energy theorem one 
finds for the total enthalpy 

DB_ldp 

Dt~ pdt [6) 

If one writes dp = a~ 2 dp+p s d s where p s denotes the derivative of the density 
with respect to the entropy, one gets from the continuity equation 

p (DB\ ds , A . 

where w denotes the mass flux. If one multiplies Crocco's vortex theorem 
([21) with the density p one obtains for the mass flux w the equation 

<9w pv D B „ 

9T-^UT + ' >VB = -' ,L + * V = K ' (5 > 

It is easy to eliminate the mass flux w from these equations. One then 
obtains an equation, which is linear in B 

„ „ „ . pv D B . d p DB ,. (9o s 

£B = V .(pVB-^— = -d,v,L +1 | r+ d 1 v 9 ,v = qM . (6) 

The operator C obtained here agrees completetely with the operator £ pot 
given in the previous section for the propagation of irrotational sound waves. 
If one inserts the sources from the equations (jlj) und ([5]), one finds 

( ° ° i ° i m „\ i , M 
= I foP'Qj: + dlv ^ sV ^ + dlv P T " j s + dlv P v x w - C 7 ) 

The sources are linear expressions in the vorticity vector and the entropy. In 
this analogy one may think of the sound as being generated from vorticity 
and entropy inhomogeneities. 

With acoustical applications in mind, the total enthalpy was first used by 
Howe [5] and recently proposed also by Doak |6J. A comparison with their 
equations shows, that the equation for B is not uniquely determined. Howe's 
convected wave operator is 

1 Dv D 1 DB 

CuoweB = AB -— ■ VB - -prr-^-prr 

a z Dt Vt a z Vt 
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and Doak's 



£r>nak B = AB 



a 



d 2 B ( d \ 

2v-+wxv+TVs-2V/i VB+vv-V-VB 



dt 2 \ dt 



The right hand sides of these equations differ from the right hand side of 
eq.(|5J. We will restrict ourselves to eq. (EJ). Notice however that the principle 
parts - i.e. those terms which contain second derivatives of B - of the 
three convected wave operators agree. This means that they agree in the 
high frequency limit of geometric acoustics and agree with the well known 
geometric acoustic theory. 

A certain simplification of eq. is possible if one requires that p and v 
fulfill the continuity equation. One then gets 

CB^.( P VB)- P §-L££. (8 ) 

The equation (jSD or (|S]) is a generalization of the wave equation. It re- 
duces to the wave equation if one assumes in flSD v = and if one assumes 
further, that p and a 2 are constant. The equation agrees then with Lighthill's 
equation in the form of Powell [10J. It seems to be a rather complicated equa- 
tion. Considering a variational principle we will however see, that equation 
(161151) is actually the simplest equation which contains a flow velocity. A sim- 
plification of the source is possible for an ideal gas. For an ideal gas p is the 
product of two functions which depend only on p and on s. Then one has 
p s = pf(s) with some function f(s) and one obtains for the sources of eq. 

m 

D ds <9v 
q tot = -divpL + p— f(s)— = -divpL-p— ■ /(s)Vs 

if one makes use of the entropy conservation law. 

The equation (H]) and (EJ) rsp. §6§ are now considered as the basic equa- 
tions of the acoustic analogy. They agree formally with the linearized equa- 
tions which describe small perturbations of a steady potential flow. The 
zeroth order equations are then given by Bo = (because of the Bernoulli 
equation) and div w = 0, i.e. the zeroth order versions of the eqs. (jl]) and 
The first order eqs. of (jl]) and §5§ are then obtained, if the fields p, a, 
and v are replaced by their zero order values. Notice that pv in eq. ([S]) then 
becomes poVo and it differs from w which becomes Wi. In that sense the 
left hand side of the eqs. (jl]) and ([5]) rsp. (Q are considered as equations 
which describe sound propagation in potential flows. This is very similar to 
Lighthill's interpretation of the wave equation as an equation which describes 
sound propagation in a medium at rest. This interpretation is valid only if 
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p, a, and v are steady fields which fulfill the equations for compressible ir- 
rotational flow. We will however not require this, as C is also well defined 
without this assumption and there are no advantages in assuming it. The 
eqs. @ and ([H]) rsp. © are exactly valid identities. They are considered in 
the following as a system of linear partial differential equations for the vari- 
ables B and w. This implies also that the right hand sides of the equations 
fll]) and §5§ rsp. (jSD are considered as the sources of the sound. They are 
related to vortices and entropy inhomogeneities. This seems reasonable if one 
wants to study the generation and propagation of sound in a potential flow. 
Then the flow consists of a superposition of an irrotational steady part and 
an unsteady part. Often one will superpose these two contributions linearly. 
We will not require that. For the steady flow one has a constant value of 
the total enthalpy B. One may then assume B to be zero. B is then solely 
related to the unsteady part and is small if this part is small. In irrotational 
flow there is a potential $. This obeys the Bernoulli equation 

$t + B = (9) 

B then differs in regions where the flow is irrotational from the temporal 
derivative of the potential only by its sign, it is however - contrary to the 
acoustic potential defined everywhere. 

This equation flS]) was originally derived in [7j. Let us now derive its 
main properties. The first important point is that C is formally self-adjoint. 
This follows from the fact that one has for arbitrary functions B and B the 
identity 

BCB - BCB = ^ + ^ (10) 

dt dx l v ' 

with 

a 2 V Dt Dt / ' 1 ^ Q x i da? J ^ 1 v > 

This equation is easy to check. It implies also that one has for a scalar 
product (f,g) defined by (f,g) = J f gd 3 xdt the relation 

(B,CB) = (B,CB) (12) 

if B and B vanish on the boundary of the integration region or decay suffi- 
ciently rapidly at infinity. 
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2.3 Reciprocity 

One may derive from the symmetry in eq. (|T2|) a reciprocity relation. To be 
specific let G(x, t,y,t') be the Green's function associated with C, i.e. 

£G(x,t,y,t') = -6(x-y)S(t-t'), with G(x, t, y, t') = for £ < t', 

(13) 

where we have assumed that G is causal, i.e. it vanishes for all times before 
the source is switched on, which occurs at t — t'. There exists also an 
advanced Green's function G adv with 

CG adv (x,t,y,t') = -5(x-y)5(t-lf), with G adv (x, t, y,t') = for t>t'. 

(14) 

It is easy to see that eq. ([6]) is invariant with respect to time reversal t — > — t 
if at the same time the sign of the velocity is reversed. Therefore time reversal 
transforms a Green's function into a Green's function. As the time reversal 
interchanges the inequalities t > t'and t > t' one has 

G(x, -t, y, -t'; -v(x, -t)) = G adv (x, t, y, t'; v(x, t)) (15) 

where we have added for clarity the function v(x, t) to the list of arguments 
of the Green's function. In addition the functions p(x, t) and a(x, t) have to 
be replaced by p(x, —t) and a(x, —t). 
One may now apply eq. fll2j) with 

£ = G(x,t,y,0 and B = G adv (x, t, z, t") 

and one obtains 

(G adv (x, t, z, t"), £G(x, i, y, 0) = (>CG adv (x, t, z, t"), G(x, t, y, 0) (16) 

Let us indicate briefly that there are no contributions from the right hand 
side of eq. ffTOj) . If the integration in eq. (TT6]) is performed over a large 
cylinder in the x,t-space which extends over a large sphere in x-space and 
over all t with T < t < Ti, one has surface contributions which are to be 
evaluated over the large sphere at t = T and at t — T\ and over the surface 
of the large sphere in x-space for all t with Tq < t < T±. If Tq is before t' and 
t" and Ti after t' and t" there are no contributions from the space integrals 
at t = Tq and at t = Ti, as at least one factor vanishes in l and in the 
/j, namely the factors containing G adv at t = T\ and those containing G at 
t = Tq. There is also no contribution from the large surface in x-space if it is 
selected so large that no signal which was generated at t — if and at x = y 
has reached this surface. Therefore eq. ffT6l) is true, and one can evaluate 
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the scalar products with the (^-functions in the eqs. (IT5|) and (fll]) easily, and 
one obtains the equation 



-G adv (y,t',z,t") = -G'(z,t",y,t') (17) 
which can with eq. (TT51) be rewritten as 

G(y, -t', z, -t"; -v(y, -t)) = G(z, t", y, t'; v(z, t)). 
This is the reciprocity principle with reversed flow. 

2.4 Variational Principle and Energy Conservation 

From the self-adjointness on can conclude the existence of a variational prin- 
ciple from which eq. fl6]) can be derived. One has 

6L = mit L = -(B,£B)-(B,q tot ) (18) 



as 



5L = ±(5B, CB) + ±(B, C5B) - (5B, q tot ) = (5B, CB - q tot ). (19) 



The Lagrangian from ([18]) can be simplified somewhat, if the second 
derivatives in L are eliminated with integration by parts. One finds then 



L 



2a 2K Dt' 2 



-(VB) 2 -q tot B dtxdt. (20) 



It seems very remarkable that this variational principle seems to be the 
simplest possible extension of the well known principle for the wave equation 
which is invariant with respect to Galilei-transformations. The Lagrangian 
in eq. (IT5|) differs from the Lagrangian of the wave equation only by the 
fact, that partial derivatives with respect to time are replaced by material 
derivatives formed with the velocity field v. 

Now it is possible to obtain from a variational principle an energy the- 
orem, rsp. an energy conservation law, if the Lagrangian density does not 
depend explicitly from the time. If I denotes the Lagrangian density from 
eq. (J20D, i.e. 

one has for the energy theorem 

d / • dl \ d ■ dl dl 

di( B dB~ ' B dB- = m (22) 
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where the time derivative on the right hand side acts only on the explicit 
time dependance in I, i.e. here in p, a, v and q to t- The variable B and its 
derivatives are to be kept constant. The time derivative on the left hand 
side of eq. (l22|) acts also on the implicit dependance in B and its derivatives 
B und B x i. Only the spatial coordinates x % are to be kept constant. An 
energy conservation law is obtained from (T2"2"j) if the Lagrangian density does 

not contain the time explicitly. In general one obtains for the energy flux 

• dl ■ dl 

Ui = B and the energy density e = B-^ — I the explicit expressions 



and 



which shows that the energy density e is positive for subsonic v. The energy 
theorem is then of the form 

oe ,. TT 

— + div U = q E n 
at 

with an energy density e and an energy source densitye q^-a- A useful relation 
is obtained if one integrates this equation over the time for finite time events 
or if one averages this equation for the case of steady phenomena. If one 
denotes the resulting quantities by an overbar, one obtains 

divU = g En . (24) 

Another important conclusion can be drawn from the energy theorem if 
one applies it to the solution of an initial value problem with vanishing right 
hand side qtot- If one assumes that the solution vanishes for large |x|, one 
obtains 

f 

— / e d 3 x = 0, 
dt J 

i.e. the toatal energy in the sound field remains constant. As it is a sum of 
positive contributions, none of these - e.g. B - can grow exponentially in 
time, i.e. instabilities cannot occur. 

The physical meaning of this energy flux becomes clearer if one considers 
an irrotational isentropic region. There on may write (1231) with the eqs. ([3]) 
and ([HD as 

U = fi(^v + pV$ t ). (25) 
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One may compare this energy flux with the flux from the Blokhintzev energy 
theorem which is valid for the propagation of irrotational sound waves in an 
irrotational mean flow in linear approximation. One has neglected quantities 
which are quadratically in the acoustical quantities. Notice that no lineariza- 
tion has been assumed in the derivation of the energy theorem (I21ll22p . It is 
insofar an exact identity, only dissipative effects have been ignored. If one 
denotes in the Blokhintzev energy flux with Ubi, the density, the speed of 
sound and the velocity of the irrotational mean flow with po, a o and vo and 
with p' and 0' the acoustic pressure and the acoustic potential, one can write 

Ubi = {-p' + v • V0')(^Ao + poW)- (26) 
Po % 

A comparison of eq. (12"6"1) with eq. (|25p shows, that both energy fluxes 
are products of two factors, where the factors of eq. (|25p are just the time 
derivatives of the factors of the Blokhintzev energy flux (1261) . If one thinks of 
the sound field as a superposition of temporal Fourier modes, one finds that 
the average energy flux consists of a superposition from fluxes of the modes. 
In the energy flux of eq. (T2"5j) all contributions contain an additional factor 
u 2 if u denotes the angular frequency of the Fourier modes. 



2.5 Solutions 

As a first application one may consider the case of constant values of the 
velocity v, density p and speed of sound a. Then vorticity and entropy 
inhomogeneities are convected with the velocity v. These inhomogeneities 
and also the total enthalpy B are then functions of x — vt only, i.e. B = B(x 
- v t), s = F s (x - v t) and u = F w (x - v t). Then eq. ()6]) leads to 

( d d d \ 

pAB = q tot , q tot = I -7^{ps-Q t + div P*v— ) + div pTV I F s + div pv x F u . 

One may then introduce B u and B s by 

P AB W = Fu und P AB S = F s 

and one obtains 

(d d d \ 

— (p s — + divp s v— ) + divpTV j B s + divpv x B w . (27) 

This shows - as one might have expected - that passively convected entropy 
and vorticity inhomogeneities do not radiate sound. In the general case these 
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quantities will not be passively convected and one needs for its determina- 
tion extra equations. For the entropy one may use the equation of entropy 
conservation 

^ = 
Dt 

and for the vorticity the Beltrami vortex theorem 

= - ■ Vv. 

Dt p p 

Here we consider especially the two-dimensional case. Then the right hand 
side of the Beltrami vortex theorem, which is related to the stretching of 
vortex lines, vanishes. 

In addition to the differential equations one needs boundary conditions. 
If one is interested in cases where the vorticity vanishes at solid walls in 
the flow region, one may use the relation ([9]) between total enthalpy B and 
potential $ and one finds that the normal component of the velocity vanishes 
at a fixed surface if the normal derivative of B vanishes there. With eq. (j23p 
one notices that the normal component of the energy flux (1231) vanishes at 
rigid walls if the normal component of the velocity v vanishes there. One 
would expect this of course. 

In addition one needs conditions of no-reflexion at the boundary of the 
computation region. We will here assume the simplest quasi-onedimensional 
condition and require there 

OB 

— = -(v + na) • VB, 

where n denotes the outer normal of the computation region. 

As a numerical example, we consider a localized vortex of radius 1 and 
of vanishing total vorticity which is convected in a flow around a circular 
cylinder of radius 1/2 and situated at x = and y = 0. The initial azimuthal 
velocity w around the center of the vortex, which is situated initially at 
x = —3, y = 0.5, is assumed to be 

w = (1 -4r 2 )(l -r 2 ) 2 , 

r denotes the distance from the vortex center. For the velocity, we assume 
an incompressible potential flow of velocity 1 in x-direction at a; = — oo. 
The density is chosen as 1, the speed of sound as 2. Initial values for B are 
obtained from eq. (127ft and are given by 

B = 0.5(y-y )(l ~ r 2 f for r < 1. 
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Figure 1. The 5-field generated by a vortex convected along a rigid 
cylinder. Grayscales correspond to \B\, white to B = 0. 

This problem is treated with the general purpose PDE-solver PDEase/2 
which runs on a PC. A grayplot of the I?-field at two different times is 
shown in figure 1. 

The left half shows a very early stage with the vortex to the left of the 
cylinder. As the vortex is inserted in an inhomogeneous velocity field, sound 
radiation begins immediately. The right half shows a later stage, where the 
vortex has approached the cylinder. As the flow Mach number is not small, 
one notices significant deviations from a dipole character. Hydrodynamic 
instabilities are not observed, but numerical small-scale errors are noticeable. 
A reliable numerical solution of eq.([6]) requires obviously more efforts. 
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